First principles calculation of uniaxial magnetic anisotropy and magnetostriction in 

strained CMR films 
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We performed first - principles relativistic full-potential lin- 
earized augmented plane wave calculations for strained tetrag- 
onal ferromagnetic La(Ba)MnC>3 with an assumed experimen- 
tal structure of thin strained tetragonal Lao.67Cao.33Mn03 
(LCMO) films grown on SrTiO 3 [001] and LaAlO 3 [001] sub- 
strates. The calculated uniaxial magnetic anisotropy energy 
(MAE) values, are in good quantitative agreement with exper- 
iment for LCMO films on SrTi03 substrate. We also analyze 
the applicability of linear magnetoelastic theory for describing 
the stain dependence of MAE, and estimate magnetostriction 
coefficient Aooi- 



I. INTRODUCTION 

Magnetic multilayers based on "colossal" magnetore- 
sistive (CMR) materials are important for many mag- 
netic applications including recording media and magne- 
toresistive sensors. The particular class of CMR materi- 
als, Lai_ x A :E Mn03 (A=Sr,Ca,Ba), is particularly inter- 
esting because they share the same basic perovskite crys- 
tal structure with many dielectrics, superconductors, and 
ferroelectrics. This stuctural similarity opens interesting 
possibilities for the epitaxial growth of CMR magnetic 
heterostructures for device applications. 

Magnetic anisotropy energy (MAE) and magnetostric- 
tion in magnetic CMR heterostructures is significantly 
different from the MAE in bulk manganites jij . Uniaxial 
magnetic anisotropy in thin Lao.67Cao.33MnC>3 (LCMO) 
films grown on SrTiO3[001] (STO) substrate has been 
measured and interpreted to be due to the strain arising 
from a film/substrate lattice mismatch Q. The possibil- 
ity of producing strained LCMO films grown on a LaA103 
(LAO) substrate with the easy magnetization axis along 
film normal has been proposed. 

In order to tackle the origin of the uniaxial MAE in 
strained CMR films we performed first - principles cal- 
culations of MAE for the strained tetragonal ferromag- 
netic La(Ba)Mn03 (L(B)MO) with an assumed experi- 
mental structure of strained LCMO-films grown on STO 
and LAO substrates. We also analyze the applicability 
of linear magnetoelastic theory for describing the stain 
dependence of MAE, and estimate magnetostriction in 
LCMO-films. 



II. COMPUTATIONAL METHOD 

The experimental data of the MAE and magnetostric- 
tion in LCMO films have been obtained from magne- 
tization curve measurements jijj]. The 58 nm thick films 
were grown on [001] oriented STO substrate using atomic 
layer by layer molecular-beam epitaxy Q . X-ray diffrac- 
tion data showed that the films have a tetragonal unit 
cell with in-plane lattice constant 7.3696 a.u. which is 
consistent with STO substrate lattice constant and per- 
pendicular lattice constant of 7.2373 a.u. This is fully 
consistent with the so-called "coherent" regime of the 
film strain in the "elastic" approximation pi : film in- 
plane lattice constant matches with substrate in-plane 
lattice constant that results in perpendicular tetragonal 
strain determined by film Poisson ratio. As a result there 
are no misfit dislocations and no thickness dependence of 
the perpendicular strain. 

An anisotropic energy density of a tetragonal ferro- 
magnetic film on non-magnetic substrate can be written 
as B: 



E/V = -K\m\ - Klm% - K^m 2 x m 2 y 
2 



(1) 



{K{mi 



where, K v terms are 2nd and 4th-order volume-type 
anisotropy constants, K s terms are surface/interface- 
type 2nd and 4th-order anisotropy constants, rn XtVtZ are 
magnetization cosines with respect to the crystal axes, 
and t is a thickness of magnetic film. For the particular 
case of LCMO films considered here t is w 150 ML and 
surface term in Eq.(l) can be safely neglected. There- 
fore, as a computational model for LCMO film grown 
on STO and LAO substrates we used strained tetragonal 
bulk La(Ba)Mn03 with the crystal structure parameters 
chosen in accordance with experiment Q. We also as- 
sume that the 4th-order terms in Eq.(l) are significantly 
smaller (as usual) than 2nd-order uniaxial anisotropy 
constant K\, and did not consider them. 

Since the measured films are ferromagnetic at low tem- 
perature for x — 0.33 Ca consentration, we considered 
the only ferromagnetic phase. We also used end-point 
ferromagnetic LaMn03 and BaMn03 instead of actual 
LCMO in spite of the fact that both LMO and BMO 
have an antiferromagnetic ground state. The use of ferro- 
magnetic LMO and BMO to approximate ferromagnetic 
LCMO films allows us to calculate the possible range of 
values for the uniaxial anisotropy in real LCMO films and 
to analyze the effect of electron concentration (changing 
Mn-site d-occupation) on MAE. 
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We used a relativistic version |6j of the full-potential 
linearized augmented plane wave method Q to obtain the 
self-consistent solutions of Kohn-Sham-Dirac equations 
and ground state charge and spin densities for the mag- 
netization directed along [001]-axis. For both strained 
LMO and BMO we used the special k-point method |7]] 
for the Brillouin zone (BZ) integrations. 348 k-points 
mesh in the l/8th irreducible part of the BZ || is em- 
ployed for the self-consistent calculations and a Gaussian 
broadening (a = 0.0014 Ry) is used for the eigenstates 
weighting. Lattice harmonics with angular momentum / 
up to 8 are used to expand the charge and spin densities 
and wavefunctions within the muffin - tin sphere. More 
than 100 plane waves per atom/spin are used as the first 
variational basis set to solve the scalar-relativistic Kohn- 
Sham equations; all occupied and empty states up to 2 
Ry above Ep are used as a second variational basis set 
to calculate spin-orbit coupling matrix elements. Self- 
consistency is achieved to within 1 x 10~ 5 e/(a.u.) 3 for 
charge and spin densities and to within 2 x 10~ 5 i?y for 
the all-electron total energy. 

The spin and orbital magnetic moments for magneti- 
zation along the [001] axis for strained LMO and BMO 
for STO and LAO substrates are shown in Table 1. The 
orbital magnetic moments are antiparallel to spin mo- 
ments, which is consistent with the atomic third Hund's 
rule for the case of less than half-filled d-shcll. There 
is a decrease of spin moment and an increase of abso- 
lute value of the orbital moment with substitution of La 
by Ba due to decrease of Mn d-occupation. Assuming 
single-site approximation for the magnetic moment we 
can estimate LCMO spin moment to be 2.85 hb and or- 
bital moment -0.011 hb for the STO substrate and 2.77 
and -0.014 fis respectively for LAO substrate. 

The MAE Q is then obtained by applying the force 
theorem to the spin - axis rotation (Tuj: from the self- 
consistent ground state charge and spin density obtained 
for the [001] spin axis, a calculation of the band structure 
for [100] spin axis orientation is performed; difference of 
the single particle eigenvalue sums is then taken to be 
the MAE. For the MAE calculations we used special k- 
points mesh which was chosen with respect to the mag- 
netic symmetry for the [100] spin axis fuJl . Since the lo- 
cal force theorem was used, the symmetry of eigenvalues 
rather than Hamiltonian was considered: for the tetrag- 
onal symmetry and magnetization directed along [100] it 
leads to exclusion of fourfold rotations with respect to 
[001] axis and leaves only eight space group symmetry 
operations (mmm space group) which are used to gener- 
ate a set of irreducible k-points in l/8th of BZ (different 
from those for self-consistent caculations). 

To achieve convergent results for the MAE, we have 
increased the number of the BZ k-points and have done 
the band calculation for both the [001] and [100] spin 
axes using the ground state charge and spin densities. 
The MAE dependence on the number of k-points for the 
LMO(STO) is shown in Fig. 1. Surprisingly good con- 
vergence (less than 1 /J,eV) was achieved for 1000 k-points 



in l/8th of BZ (8000 in full BZ). This set of k-points was 
then used for all MAE calculations. 



III. RESULTS AND DISCUSSION 

The calculated uniaxial MAE for all four cases 
(L(B)MO(STO) and L(B)MO(LAO)) are shown in Table 
II. The calculated MAE values are in very good agree- 
ment with experiment for LCMO films (—56 ^eV) on a 
STO substrate. This is direct numerical evidence of the 
magnetoelastic origin of the uniaxial MAE in strained 
LCMO films and supports quantitatively an interpreta- 
tion of experiment given in Ref. §. For LMO (LAO) and 
BMO (LAO) we calculated the MAE to be positive. It 
is in qualitative agreement with very recent experiments 
for Lao.8Cao.2Mn03 films grown on LAO substrate Q. 

In addition to calculated magnetocrystalline MAE 
caused by spin-orbit interaction we have to take into ac- 
count demagnetization energy (DE) which is due to "clas- 
sic" magnetic dipole-dipole interaction. For the ferro- 
magnetic film the DE can be evaluated numerically [ fLs] . 
We calculated the DE using the results of Ref. Q and 
taking into account that magnetic moment in L(B)MO 
is located at Mn-site (cf. Table I.). 

The values of DE for all four cases (L(B)MO(STO) 
and L(B)MO(LAO)) are shown in Table II. (We took 
into account the only in-plane dipole-dipole interaction 
between Mn-sites.) For the case of L(B)MO(STO) both 
MAE and DE keep magnetization in the film plane. For 
the case of LMO (STO) the sum of MAE and DE is neg- 
ative (—28 neV) and magnetization is in the film plane. 
For the BMO(STO) the total uniaxial MAE is positive 
(42 /J,eV) and magnetization is directed along the film 
normal. 

Assuming single-site approximation for the magnetic 
anisotropy we can estimate total uniaxial MAE ss 
—7 iieV per cell for LCMO film on LAO substrate. 
It means that the absolute value of the magnetoelastic 
MAE is very close to the absolute value of the demagne- 
tization field for LCMO film on LAO substrate. 

To check the applicability of linear magnetoelastic the- 
ory we map our results onto the usual anisotropic mag- 
netoelastic energy dependence on strain: 

E = B\ (e z — eo) m 2 z + const , (2) 

where eo is a biaxial (in-plane) strain due to the 
film/substrate lattice mismatch, e z is a uniaxial (along 
[001]) strain, B\ is a magnetoelastic coefficient, and m z is 
a magnetization cosine with respect to the z([001])-axis. 

Using Eq.(||), the values of strain Q and MAE, we 
calculate the values of B x for L(B)MO for STO and 
LAO substrates (cf. Tab. II). There is quite pro- 
nounced variation of B\ for LMO for the different sub- 
strates. It means that for LMO there is significant devi- 
ation from linear MAE dependence on strain (Eq.(||)), 
and linear theory is rather qualitative. For the case 
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of BMO the variation of B\ is smaller, therefore lin- 
ear theory is more reliable. Note, that calculated val- 
ues of magnetoelasic coefficient B\ are in very reason- 
able agreement with the experimentally derived value 
-6.7 x 10 7 erg/cm 3 (-2.44 meV per cell) for LCMO 
films on STO substrate. 

In spite of the fact that linear theory gives the only 
a semi-quantitative description of uniaxial MAE depen- 
dence on strain, our calculations show that one can still 
use it for estimation of the sign and order of magnitude of 
the magnetostriction coefficient Aooi- To determine Aooi 
one has to minimize the sum of Eq. (0) and elastic energy 

1 2 

E e lastic = 2 Clie z + 2c 12 e z e + Const , (3) 

where, en and C12 are elastic moduli, with respect to 
e 2 , with eo fixed by substrate. It leads to the following 
expression for the magnetostriction constant Aooi (see for 
details Ref. §): 

Aooi = - \ — (4) 
3 en 

Assuming c\\ = 2 x I0 12 erg/cm 3 yj] and using cal- 
culated values of B\ (cf. Tab. II) we calculate Aooi for 
L(B)MO films on STO and LAO substrates (cf. Table 
III). Our calculated values of Aooi agree in sign and or- 
der of magnitude with the experimentally derived value 
of 7 x 10~ 5 . However, we have to note that it is unclear 
what definition for Aooi in terms of B\ and en was used 
in Ref. [Q : using the data of Ref. |jl| and Eq. (||) we cal- 
culated A oi = 2.23 x I0 -5 in very good agreement with 
our results. 



IV. CONCLUSION 

To summarize, we have shown quantitatively that 
the observed uniaxial anisotropy in LCMO films grown 
on STO substrate is caused by strain arising from 
film/substrate lattice constant mismatch. It is also 
shown that the magnetoelastic MAE is positive for 
LCMO films on LAO substrate and absolute value of 
MAE is close to the absolute value of demagnetization en- 
ergy. By varying the strain arising from a film/substrate 
lattice constant mismatch one can control the uniaxial 
anizotropy for LCMO films in order to overcome the 
demagnetization energy and provide [001] spontaneous 
magnetization direction. 

I am grateful to W. E. Pickett, A. J. Freeman and O. 
N. Mryasov for helpful discussion. This research was sup- 
ported by the National Science foundation Grant DMR- 
9802076. 
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TABLE I. Experimental lattice constant (a, in a.u.) and 
c/a ratio for tetragonal strained LCMO films on STO and 
LAO substrates; spin (M s ) and orbital (Mi) magnetic mo- 
ments for LMO and BMO for the [001] spin direction. 



Substrate: 


STO 




LAO 




a — D 


< .ooy / 




7.1618 




c/a 


0.982 




1.047 




Film: 


LMO 


BMO 


LMO 


BMO 


M s (Mn) 


2.996 


2.563 


2.880 


2.542 


M s (Total) 


3.145 


2.727 


3.022 


2.712 


Mi (Mn) 


-0.004 


-0.023 


-0.007 


-0.023 


Mi (Total) 


-0.002 


-0.021 


-0.006 


-0.024 


TABLE II. 


Biaxial (eo 


and uniaxial (e z ) strains, uniaxial 


MAE (fxeV), 


demagnetization energy 


(DE) (peV) 


Bi (meV) 


for LMO and BMO films 


on STO anc 


I LAO substrates. 


Substrate: 


STO 




LAO 




eo 


0.008 




-0.021 




e z 


-0.0101 




0.0258 




Film: 


LMO 


BMO 


LMO 


BMO 


MAE 


-40.9 


-53.3 


37.7 


100.2 


DE 


-65.7 


-48.7 


-66.1 


-58.3 


Bi 


-2.259 


-2.944 


-0.806 


-2.141 



TABLE III. Aooi x 10 5 for LMO and BMO films on STO 



and LAO substrates. 



Substrate: 


STO 




LAO 




Film: 


LMO 


BMO 


LMO 


BMO 


Aooi 


2.1 


2.74 


0.76 


2.01 



50 




500 1000 1500 2000 

k-points, 1/8 BZ 

FIG. 1. The calculated uniaxial MAE dependence on 
k-points number in l/8th irreducible part of the BZ. 
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